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1. Introduction 

It is an old idea that \wg&-N QCD might be exactly reformulated as a string theory. Although 
this relation has never been made precise, the results obtained in lattice gauge theory leave little 
doubt that there is a connection between these two theories. One of these results that leads to such 
a conclusion is the energy spectrum of the flux-tube. A flux-tube in SU {N) gauge theories with 
length / much larger than its width is expected to be described by an effective low energy string 
theory. This can be verified by studying the energy spectrum of the flux-tube and comparing it to 
the effective string theory prediction. 

In this work we focus on the spectrum of the closed fundamental flux-tube in pure SU {N) 
gauge theories in D = 2 + 1 dimensions. We calculate the energies of the lowest lying ~30 states. 
This allows to extend the comparison with theoretical predictions to states with more quantum 
numbers and a nontrivial degeneracy structure. More precisely, we want to check whether the 
closed flux-tube can be described by the Nambu-Goto (NG) string model |[T]l, and, if so, how good 
a description it is. 

The flux-tubes that we study have lengths that range from ~ 0.65 fm to ~ 2.60fm and the 
lattices we use have spacings a ~ 0.04 fm and a ~ 0.08 fm (Despite working in pure gauge theories 
and in D = 2 + 1, we define Ifm through the convention a = (440MeV)^). The gauge groups that 
we study have N = 3,6. The use ofN>3 will suppress any mixing amplitudes (eg. glueball-string 
mixing) that cannot be described by a simple low-energy effective string theory model. 

The study of confining flux-tubes with lattice techniques has been an active field for the past 
three decades, and we refer the reader to some recent papers For more details on the calculation 
and relevant references see our longer write-up 

2. Lattice construction 

We define our gauge theory on a three-dimensional periodic Euclidean space-time lattice with 
L X L± X Lj sites. It is important to mention that in order to minimize the finite volume effects, we 
increase Lj_ and Lj as we decrease the length of the flux-tube. For the calculation of the physical 
observables we perform Monte-Carlo simulations using the standard Wilson plaquette action: 



The bare coupling j3 is related to the dimensionful coupling through limQ^o jS = 2N/ag^. In the 
laige-N limit, the 't Hooft coupling A = g^N is kept fixed, and so we must scale j3 = 2N^/X oc 
in order to keep the lattice spacing fixed. The simulation we use mixes standard heat-bath and 
over-relaxation steps in the ratio 1 : 4. These are implemented by updating SU (2) subgroups using 
the Cabibbo-Marinari algorithm. 

We have calculated the string spectrum for the case of SU{3) with a ~ 0.04,0.08 fm and 
SU (6) with a ~ 0.08 fm. For the case of SU (3) and a ~ 0.08 fm, the string lengths ranged between 
~ 0.65fm and ~ 2.60fm and for the two other cases between ~ 0.95fm and ~ 2.0fm. 

3. General strategy 

We calculate the energies of flux-tubes that are closed around a spatial torus using the correla- 
tors of suitably smeared Polyakov loops that wind once around the corresponding spatial tori and 
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that have vanishing transverse momentum. This is a standard technique with smearing/blocking de- 
signed to enhance the projection of our operators onto the physical states. We classify our operators 
using the quantum numbers of transverse parity in P = it and winding momentum ^ = 0, ±1 , ±2, . . . 
in units of 2n/l. 

For each combination of these quantum numbers we construct the full correlation matrix of 
operators and use it to obtain best estimates for the string states using a variational method applied 
to the transfer matrix f = e^"^ . 

Since we are mostly interested in the excited states, it makes sense to introduce transverse 
deformations in the simple Polyakov loop. Using this procedure we can increase the number of 
different operators to 80-200 and construct Polyakov loops described by the quantum numbers of 
P and q. We present the different paths used in the calculation in Table [T[ 
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Table 1: The lattice paths used in the construction of Polyakov loops in this work. Our set of operators 
can be divided into three subsets: a. The simple line operator (1) in several smearing/blocking levels, b. 
The wave-Uke operator (2) whose number depends upon L, L±, and the smearing/blocking level, c. The 
pulse-like operators (3-15) in several different smearing/blocking levels. 



4. Theoretical expectations: The spectrum of the Nambu-Goto String model 

The action of the Nambu-Goto (NG) model ^ is the area of the worldsheet swept by the prop- 
agation of the string. This model is quantum-mechanically inconsistent due to the Weyl anomaly 
(D 7^ 26), but since this anomaly is suppressed for long strings |Q] it can still be considered as an 
effective low energy model. The classical ground state of the model is the worldsheet configuration 
with the minimal area, and the transverse area fluctuation around this ground state constitute the 
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quantum spectrum of the model. These fluctuations are usually referred to as left/right movers, 
depending on the momentum that they carry in the direction parallel to the string axis. 

The single string states can be characterised by the winding number w which indicates how 
many times the string winds around the torus, by the occupation number and the energy k 

of the left and right movers and also by the center of mass momentum Pc.m.- By projecting to zero 
transverse momentum we are left only with the momentum along the string axis which is quantized 

in units of 2nq/l with q = 0, ±1, ±2, These quanta are not independent of ni R and obey the 

level matching constraint Nl — Nr = qw, where A^^(s) enumerates the momentum contribution of 
the left(right) movers in a certain state as follows: 

Nl=Z E Mk)k and iV« = £ £ nR{k')k' (4.1) 

k>OiiL{k)>0 k'>OnR{k')>0 

The string states can be characterised as irreducible representations of the SO{D — 2) group that 
rotates the spatial directions transverse to the string axis. In our D = 2 + 1 case this group becomes 
the transverse parity with eigenvalues P = (— l)£'=i"i('^')+^j=i"'^('^j). Finally, the energy of a closed- 
string state described by the above quantum numbers for any D is given by the following relation: 

^i*,„, - C^'^f + 8.. - ^^^) + (^) ' . (4.2) 

5. Theoretical expectations: Effective string theories 

Since in 2+1 dimensions the NG string is at best an effective low-energy string theory, it makes 
sense to generalise it and write the most general form of an effective string action consistent with 
the symmetries of the flux-tube system. This was first done in the early eighties for the case of 
w = 1 and ^ = in 1^, and the spectrum obtained for any D was given by: 

£. = aZ + ^(n-^)+^(l/z2), (5.1) 

where n = 0, 1,2, .... The second term in the above formula is known as the Liischer term and 
is expected to be universal. One can easily show that the NG model obeys this universality by 



expanding the square-root of Eq. (^) to leading order in I /I. 

Recently, the work [||] was extended in [^, where the authors used an open-closed string 
duality of the effective string theory. Using this duality they showed that for any D the is 
absent from Eq. (^?T|), and that in D = 2 + 1 the has a universal coefficient. Consequently, 



in 2 + 1 dimensions Eq. (5. 1 ) is extended to: 



47r / 1 \ SttV 1 \ 2 



£. = a/+-(«--l-^(.--) +^{l/n (5.2) 



or in the equivalent form: 



E^, = {alf + S7ia{n-^]+^{l/l'). (5.3) 



The form Eq. (5.3) is particularly convenient since the two first terms on the right hand side are the 



predictions of the NG model that we find to be a very good approximation. 
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Motivated by this recent development, we decided to fit our data using tlie following ansatz: 



p ' 



P>3 



(5.4) 



where e1;,q is the Nambu Goto prediction given by Eq. (12) for w = l,n = Nl = Nr and where C. 



are dimensionless coefficients that in general can depend on the quantum numbers of the state. 



6. Results 

In this section we present our results from the calculation. The results were obtained for 
SU{3) with p = 21.00,40.00 and for SU{6) with j3 = 90.00. We have focused on strings with 
1 .4 < l^/a < 5.5. All energies that we present were obtained from single cosh fits to the correlation 
functions of our 'best' operators - see discussion in Section |[ 

We present the results in Figs. [l|,Q The lines are the predictions of the NG model. The 
string tensions used for these predictions were extracted from the ground state energies of the 
^ = calculations with the use of the ansatz Eq. (5.4) for p = 3. The fitting parameters are the 



dimenionless C3 and the lattice spacing in physical units a^/a. We present these in Table ^ 



Gauge group 


P 




C3 


Confidence level 


SU{3) 


j3 =21.00 
j3 = 40.00 


0.030258(26) 
0.007577(13) 


0.160(21) 
0.05(31) 


69% 
15% 


SU{6) 


p = 90.00 


0.029559(36) 


0.04(21) 


88% 



Table 2: The parameters a a and C3 in the fit. 



According to our results the NG predictions are very good approximations to the flux-tube 



spectrum and deviate from our data only at the level of ~ 2% when l^/a ^ 4.2. In contrast to that, a 



comparison of our data to Eq. (5.1 ) and/or Eq.( 5.2) fails even when l^/a ^ 4. This is an important 



point that tells us that the higher 0{l/l^) terms in Eq. ( |5.2[ ) are significant for the excited states at 
these lengths, and that they are captured quite well by Eq. ( pT2| ) (to the level of ~ 2%, a deviation 
that may or may not be due to a percent level systematic errors that we did not aim to control for the 
excited states). Finally, it is important to remark that the degeneracy pattern predicted by the NG 
model is seen from our data. For example, the second energy level is fourfold degenerate at large-/. 
This degeneracy includes two positive parity states and two negative parity states, and these start 
sphtting significantly once l^/a < 3. 

Next, we performed fits of the energy of the excited states. In the case of the first excited 



energy level, where there is only one state per level, we used the fitting ansatz Eq. ( p.4[ ). In the 
case of the second excited energy level, where for each parity there are two states, we fitted the 
difference between the energies squared of these states. The fits showed that to unambiguously 
determine the power p, and test the Liischer-Weisz prediction of [^, we need statistical errors 
which are at least 2-3 times smaller than the ones our data has, and a simultaneous control of any 
systematic errors that may be important at the level of a few percents accuracy. 

Finally, we move to the non-zero longitudinal momentum q y^O calculation. In Section 4 we 
have mentioned that the number of left and right movers in the NG prediction is constrained by the 
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Figure 1: The energies of the first three energy levels divided by ^/o as a function of l^/o. The three lines 
are the NG predictions for the ground state and the excited states. The green lines are the NG predictions 
expanded to leading order Eq. (5.1) and next-to-leading order Eq. (5.2) for n ~ I. Left panel: We present 
the results for the case of SU{3) for two different values of j3 (two different lattice spacings). For j3 = 21, 
positive parity states are presented in red and negative parity states in blue. For )3 — 40, positive parity states 
are presented in cyan and negative parity states in pink. Right panel: We present the results for the case of 
SU(6) with j3 = 90. Positive parity states are presented in red and negative parity states in blue. 



level matching condition Nl — Nr = qw. The comparison of this prediction to our data is presented 

in Fig. ^ where we plot \J E'^ /o — {l%ql ^/olf' as a function of l^fd. It is clearly seen that our 
data is very well described by Eq. ( [4.2[ ). 



7. Summary 

We have calculated the energy spectrum of closed strings in the fundamental representation 
of SU{N) gauge theories in 2-1-1 dimensions. To perform this calculation we constructed a basis 
of ~ 80 — 200 operators for each configuration of quantum numbers. The use of this large ba- 
sis of operators is convenient for two reasons. Firstly, it is possible to extract masses of excited 
states with high confidence since it increases the overlaps of our lattice operators onto the physical 
states (compared to what we observe using the simple line operator in 5 smearing/blocking levels). 
Secondly, it enables us to study states with quantum numbers like the transverse parity P and the 
longitudinal momentum q. 

Comparing our results to different theoretical predictions, we find that the agreement with the 



NG prediction in Eq. (4.2) is very good, including the expected degeneracy pattern. This agreement 



is in striking contrast to what we find when we simply compare to the Luscher term as in Eq. ( p.l[ ) 
or to the Luscher- Weisz prediction in Eq. ( |5.2[ ). While for the ground states these describe our 
data already at l^/G ~ 3 where they are indistinguishable from the full NG formula, for the excited 
states the situation is completely different. In particular, whereas the NG prediction works well for 



the first excited state already at l^/a ~ 3, the predictions of Eq. (5.1) and Eq. (]5.2D are still very far 
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Figure 2: y E^/ C ^ {^TZq/y/al) as a function of ly/a for the ground state and excited states of the non- 
zero momentum along the string direction. The four Hnes present the NG predictions Eq. (4.2). For Nr — 
3,Nl = i and q = 2, the expected degeneracy is three, which agrees with our excited state data, i.e. two(one) 
states with positive(negative) parity in red(blue). 

from the data (see Figure. |I|). This demonstrates that a confining flux-tube can be described by a 
covariant string, with small or moderate corrections down to very short lengths. 
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